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ABSTRACT 


A set of diffuse interstellar clouds in the inner Galaxy within a few hundred 
pc of the Galactic plane has been observed at an angular resolution of ~ TO 
combining data from the NRAO Green Bank Telescope and the Very Large Array. 
At the distance of the clouds the linear resolution ranges from ~ 1.9 pc to ~ 
2.8 pc. These clouds have been selected to be somewhat out of the Galactic plane 
and are thus not confused with unrelated emission, but in other respects they 
are a Galactic population. They are located near the tangent points in the inner 
Galaxy, and thus at a quantifiable distance: 2.3 < R < 6.0 kpc from the Galactic 
Center, and —1000 < z < +610 pc from the Galactic plane. These are the 
first images of the diffuse neutral H I clouds that may constitute a considerable 
fraction of the ISM. Peak H I column densities range from Ahi= 0.8 — 2.9 x 
10 20 cm -2 . Cloud diameters vary between about 10 and 100 pc, and their 
HI mass spans the range from less than a hundred to a few thousands M 0 . The 
clouds show no morphological consistency of any kind except that their shapes 
are highly irregular. One cloud may lie within the hot wind from the nucleus 
of the Galaxy, and some clouds show evidence of two distinct thermal phases as 
would be expected from equilibrium models of the interstellar medium. 

Subject headings: ISM: structure - ISM: clouds - ISM: atoms - ISM: general - Galaxy: 
disk - radio lines: ISM 
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1. Introduction 


The concept of a diffuse interstellar cloud is more than 50 years old, yet there are few 
observations that support the most basic aspects of the standard picture. The strongest 
evidence for discrete clouds is kinematic: there are usually distinct absorption lines at 


different velocities in spectra toward stars (e.g., Munch (1952); Hobbs (1978); Redfield & 
Linsky (2008)). But spectral features can be produced not only by spatial structures, i.e., 


clouds, but in a continuous turbulent medium as well (Lazarian & Pogosyan 2000). In 


contrast to the often well-defined clouds of molecular gas like the Infrared-dark Clouds (e.g., 


Rathborne et al. 2010), there is little support for the existence of discrete clouds in 21cm 


emission observations, which suggest instead that the atomic interstellar medium (ISM) 
consists of fragments of filaments and “blobby sheets”, many of which may be a consequence 


of turbulence (Kulkarni & Heiles 1987; Dickey & Lockman 1990; Miville-Deschenes et al. 


2003; Heiles & Troland 2003; Kalberla & Kerp 2009). In most direction 21cm H I emission 


maps are highly confused, leading to considerable ambiguity in determining the morphology 
and boundries of interstellar clouds, even assuming that they do exist. A mammoth study 
of Na I and Ca II absorption lines toward nearly 2000 stars within 800 pc of the Sun by 


Lallement et al. (2003) and Welsh et al. (2010) has revealed the three-dimensional structure 


of the local interstellar medium, “cell-like cavity structures” a fragmented “wall” of neutral 
gas, and what appears to be clouds “physically linked to the wall of denser gas”, but it is 
difficult to know how to generalize this result to the broader ISM. 

The situation is quite different, however, in the lower halo of the inner Galaxy, where 
there is a population of discreet H I clouds whose velocities are consistent with circular 
rotation, but whose location several hundred pc from the place separates them from 


unrelated emission (Lockman 2002, 2004) b Similar clouds can be seen at low Galactic 


1 First detections of a few prominent representatives of this population date back to 


Prata 
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latitude when their random velocity is large enough to remove confusion (Stil et al. 2006); 


others are detected in the outer Galaxy (Strasser et al.||2007; Stanimirovic et ah|2006; Dedes 


& Kalberla 2010). The clouds in the inner Galaxy are likely the product of H I supershells 


as their abundance and scale height are linked to the large-scale pattern of star formation 


in the disk (Ford et al. 2008, 2010). 


While many aspects of these “disk-halo” clouds are poorly understood, they can be 
used as test particles sensitive to the physical conditions in their surroundings, and thus 
give information about interstellar processes not easily gotten from the highly blended 
spectra typical of most observations. The disk-halo clouds in the inner Galaxy are so 
abundant that a number of them he near the terminal velocity in their direction, and thus 
near the tangent point, whose distance is determined from simple geometry. Their location, 
mass and size can be estimated with quantifiable errors. 


There have been numerous theoretical studies of the expected properties of the diffuse 
ISM as a function of location in the Galaxy, distance from the Galactic plane, sources of 
heating, etc. Strong theoretical predictions have been made, especially about the existence 
of two thermal phases in pressure equilibrium (Field et al. 1969; Wolfire et al. 1995a|b; 


Koyama & Ostriker 2009). The disk-halo clouds offer the perfect laboratories for testing 


these predictions. 


We selected a set of disk-halo clouds using observations with the Robert C. Byrd Green 
Bank Telescope (GBT) and measured them with the Very Large Array (VLA) in three 
different array configurations. The clouds were selected to cover a range of longitude and 
latitude and to be located near the tangent points of the inner Galaxy, and thus at a known 


(1964), Simonson (1971), and Lockman (1984). A very detailed history and bibliography of 
both observational and theoretical early studies of interstellar clouds and H I halo can be 


found in Chapter 1 of Pidopryhora (2006). 
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distance. 


This paper is the first in a series about these clouds. Here we discuss the observations 
and data reduction for the GBT and VLA D-array data only, taken at an angular resolution 
of 019-115 providing a linear resolution of 1.9 to 2.8 pc. We concentrate on understanding 
all sources of uncertainty. 


2. Selection of Targets 


Targeted 21cm H I surveys of regions in the inner Galaxy made with the GBT provided 


a list of diffuse clouds that might be suitable for high-resolution imaging (Lockman 2002 


2004). From these we selected a set using the following criteria: 1) The clouds have 
LSR velocities at or beyond the terminal velocity in their direction ensuring that their 
distance could be determined (see jj |2.1| ); 2) the clouds cover a range of Galactic longitude 
and latitude ensuring that different environments were probed, though all are in the first 
quadrant of Galactic longitude; 3) the clouds are relatively isolated in position and velocity 
to minimize potential confusion; 4) their 21cm emission as observed with the GBT is bright 
enough to be detectable with the VLA in a few hours. An example of the GBT observations 
used to select the clouds is given in Figure [l] Table [l] gives the cloud designation, field 
centers, and the 21cm line peak brightness temperature, FWHM, and velocity as determined 
from the GBT observations. 


2.1. Distance 


The kinematics of the disk-halo cloud population studied here is dominated 
by the circular rotation of the Milky Way with a cloud-cloud velocity dispersion 
a cc ~ 16 km s -1 (Lockman 2002; Ford et al. 2008, 2010). Toward the inner Galaxy, the 
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maximum velocity permitted by Galactic rotation at b = 0° is called the terminal velocity, 
V, and arises from the tangent point where the distance from the Galactic center is 


R t = Rq sin L Here we use Rq = 8.5 kpc, the IAU recommended value (Kerr & Lynden-Bell 


1986). The terminal velocity can be measured from observations of species such as H I or 


CO or can be approximated with a rotation curve (e.g., |Burton fe Liszt (1993); Clemens 


(1985); McClure-Griffiths & Dickey (2007); Dickey| ( |2013| ) ) . In the first quadrant of 
Galactic longitude V t > 0 and an object with Vlsr > V must thus he near the tangent 


point where its distance is known from geometry. In the current sample all clouds have 
velocities Vlsr V V so we calculate the tangent-point distance projected on the Galactic 
plane d p = R 0 cos t. A cloud’s distance from the Galactic plane is then z = d p tan b and the 

l 

distance to the cloud center d = [d? v + z 2 ) 2. 


2.1.1. Distance Uncertainties 


An estimate of the uncertainties in a cloud’s distance can be derived from the change 
in distance that would correspond to a change in the cloud’s Vlsr of a cc in that particular 
direction for an assumed rotation curve. As almost all of the clouds have Vlsr > V, in 
some cases by as much as 2cr cc , we conservatively estimate errors by the change in distance 
were the cloud to have a velocity (V — u cc ) or (Vlsr — cr cc ), whichever produces the larger 
change in distance. Errors calculated this way using the rotation curves of Burton & Liszt| 
(1993) and Dickey (2013) agree to within a few percent, and are given in Table [lj We note 
that the greater a cloud’s velocity beyond V t , the more likely that it lies near the tangent 


point (Ford et al. 2010). Thus for many of our clouds, especially G44.8 — 7.0 which lies 


nearly 30 km s 1 past V, our error estimates are probably overstated. 

The final two columns of Table [l] give derived distances and errors as well as the 
distance of each cloud from the Galactic plane. 
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3. Observations and Basic Data Reduction 


3.1. Green Bank Telescope Observations 


The GBT was used to map H I emission around the disk-halo clouds to measure their 
overall properties, to determine whether they might be suitable for high-resolution imaging, 


and to provide the short spacing data for image reconstruction (see e. g., Stanimirovic et al. 


(1999)). Maps were made over an area around each cloud of 1°5 x 1?5 or 2° x 2° depending 
on the extent of the cloud. Spectra were taken every 3(5 in Galactic longitude and latitude, 
somewhat finer than the Nyquist sampling interval for the GBT’s 9(1 beam (FWHM). 
Observations were repeated several times over a period of a few months. In-band frequency 
switching gave a useable velocity coverage of 400 km s -1 around zero velocity (LSR) at a 
channel spacing of 0.16 km s -1 . 


3.2. Very Large Array Observations 


A sample of 20 H I disk-halo clouds was observed in 21 cm line emission spectroscopy 
with the Very Large Array (VLA) in D configuration during 2003 and 2004. The spectra 
had 256 frequency channels separated by 0.64 km s _1 in equivalent Doppler velocity centered 


on the peak velocity of each cloud (Pidopryhora et al. 2004; Lockman & Pidopryhora 


2005; Pidopryhora 2006). Fifteen of the clouds were observed in single pointings, and 


5 as mosaics. Two of these were also observed in C configuration with the identical 
spectroscopic setup (Pidopryhora et al. (2012)) and two more (G21.2 + 2.2 and G35.6 + 3.9) 
were observed in B configuration at both H I and OH frequencies in an attempt to detect 
absorption against bright background continuum sources. Here we present results on just 
the ten clouds observed only in D configuration with single pointings, deferring discussion 
of the others to a later paper. Parameters of the VLA observations are given in Table [2| 
An example of the un-coverage for one of the clouds is given in Fig. [2| 
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3.3. Green Bank Telescope Data Reduction 


The spectra were calibrated and corrected for stray radiation as described in Boothroyd 


et al. (2011), and a second-order polynomial was fit to emission-free regions of each spectrum 


to correct for residual instrumental effects. For each cloud the data were assembled into 
a cube on a 105" grid. There was occasional narrow-band interference that was stable in 
frequency, so spectra were interpolated over the affected channels. The final GBT data 
cubes had a brightness temperature noise ~ 0.1 K in a 0.16 km s -1 channel. 

Each GBT image cube was converted to the same coordinate system as the VLA, 
cropped to fit the exact VLA field size and interpolated to the matching grid and sequence 


of spectral channels with the Miriad (Sault et al. 1995) task REGRID. The GBT data were 


taken while the telescope was moving and thus have an effective resolution of ~ 916 x 911, 
with the major axis along Galactic longitude, the scanning direction. In the VLA’s 
equatorial coordinate system this resulted in slightly different beam position angles for each 
field. In order to smooth out gridding artifacts each GBT image was also convolved with 
a circular beam function of 200", approximately 1/3 of the original beam size, so the final 
GBT angular resolution is approximately 10' FWHM. 


3.4. Very Large Array Data Reduction 

After calibration and study of preliminary dirty images of each VLA field, the 
continuum was subtracted in the uv domain using AIPS UVLIN task based on selected 
line-free channels. Naturally-weighted dirty images of continuum-free data were then 
cleaned channel by channel with AIPS SDCLN with no cleaning mask applied. The residual 
flux threshold was set to 0.7 mJy/beam for all clouds, corresponding to 0.2-0.3 a. In the 
final step the correction for the VLA primary beam was applied to the clean image cubes 
with the AIPS task PBCOR. The imaging synthesized beam size was different for each field 











as described below; their values are given in Table [3| 


4. Further Data Processing, Noise Levels and Errors 
4.1. Combining the Interferometric and Single-Dish Data 


For several reasons we have chosen to use Miriad’s IMMERGE to combine the 
interferometric and single-dish data: 1) it uses a well understood algorithm that is easy to 
control; 2) it runs quickly and is easily applied to large image cubes; 3) we have developed 
a calibration technique described below that ensures accuracy of the results. 


Another approach was also tried for two clouds not of the set described in this paper 
(Pidopryhora et al. 2012): using a maximum entropy (MEM) algorithm (AIPS VTESS 
task) with the GBT image used as the default, but it was found to be less efficient. For a 


detailed review of all possible methods see Stanimirovic (1999). 


Miriad’s task IMMERGE uses a linear method sometimes known as ‘feathering’ to 


combine interferometric and single-dish data (Sault & Killeen 2008; Stanimirovic 1999) 


Essentially this is just merging clean interferometric and single-dish images after Fourier- 
transforming them into the uv domain, the result covering the whole combined spatial 
frequency range. For this procedure to be meaningful the Fourier images of interferometric 
and single-dish data should match within the overlapping ranges of their common spatial 
frequencies. Due to the completely different natures of the two original datasets and thus 
unavoidable discrepancies, this has to be ensured by varying their common calibration 
scaling factor / ca i , which is the main control parameter of the IMMERGE task. In the 
case that the calibration of both data sets was done properly this factor should be close 
to unity, but its exact value has to be determined empirically in each particular case. If a 
compact source of 21cm emission were present in the field of view, both unresolved by the 
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VLA and unconfused with other emission by the GBT, determining f ca \ would be as simple 
as dividing the interferometric flux density of this source by its single-dish flux density. 
Unfortunately, such sources are rare and not present in our data, so a more complex 
strategy of determining / ca j was used. 


There is another important control parameter called ‘tapering’. If IMMERGE tapering 
is applied, the Fourier image of the interferometric data is smoothly continued into the 


low spatial-frequency region (Sault & Killeen 2008). This fixes possible edge effects in the 
Fourier transformation, but introduces an additional non-linear distortion of the data. As 
well as determining the best value of f ca \ , it is necessary to decide if tapering should be 
applied. 


4.2. Derivation of the Optimal f ca \ 

IMMERGE has a built-in method of matching the two datasets and deriving f ca \ , 
provided that the overlapping range of spatial frequencies is defined. Figure [3] shows 
examples of application of this method. The interferometric image is convolved with the 
single-dish beam, then both images are Fourier-transformed and compared to each other 
within the overlapping spatial frequency range. The value of f ca \ is selected by scaling the 
single-dish data until the slope of the line fit to the data is exactly 1. Based on several 
trials varying IMMERGE parameters with different samples taken from our data, we have 
determined that: 1) tapering distorts the data and makes a reasonable linear fit impossible; 
2) due to a large scatter of the values the best linear fit usually does not characterize the 
data well and f ca \ derived from it should be treated only as an estimate. Thus for our 
purposes we use IMMERGE without tapering and we determine the optimal f ca \ by other 
methods. 

It should be noted that describing the discrepancy between single-dish and 
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interferometric data with a single linear parameter is only an approximation of very 
complex behavior. One immediately notices that the optimal / ca j seems to be different for 


different channels (Stanimirovic 1999). In particular, it is sensitive to the signal-to-noise 
ratio and spatial signal distribution of brightness in each channel. This is understandable 
as ideally / ca i should be determined from a point source unresolved by both telescopes. For 
a cloud comprised of diffuse gas that smoothly fills the held of view and is detectable by the 
single-dish but is completely invisible to the interferometer, the estimated / ca i approaches 
zero. In the cloud data we usually have some combination of these two extreme cases, 
and what seems to be the optimal / ca i may fluctuate significantly. But allowing it to vary 
with frequency would introduce an unknown non-linear distortion. Based on the general 
assumptions of this model, and as all clouds observed only with the VLA D-array present 
the same variety of interferometric data (similar beam sizes, noise levels, m>-coverage 
etc.), we sought a single value of / ca i that would work well not only for every channel of a 
particular cloud, but for all ten clouds. 


One criterion for testing the goodness of a particular value of / ca i is that the spectrum 
of the merged cube averaged over an area somewhat larger than the GBT beam, but small 
enough not to be effected by the VLA primary beam pattern, should match the average 
taken over the same region in the single-dish data alone. We have found that averaging 
over an area 1513 x 1513 at the center of the held works well for all clouds observed with 
the VLA D-array. Requiring that the average spectrum over this area given by IMMERGE 
matches the single-dish mean profile is sufficient to derive an optimal / ca i = 0.87. In fact, 
this method may be preferable to any other since it directly preserves the single-dish flux. 

Figure [4] shows the comparison of mean line profiles of the VLA, the GBT and the best 
IMMERGE combined images for all 10 targets using the single value of / ca i = 0.87, which 
was adopted for all subsequent work. 
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4.3. Testing the Accuracy of the Image Produced by IMMERGE 


The average difference between the observed GBT and final IMMERGE line profiles 
can be calculated 


(Pgbt(*) — Pimmerge(*)) 

i 

J2 -Pgbt(^) 

% 

where Pgbt(*) and Pimmerge(z) are the spectral values for individual channels and the 
summing is done over velocities of interest for each cloud. Figure [5] shows values of (D) 
for all 10 clouds. Taking into account the peculiarities of profiles that make their exact 
match impossible, e. g., the presence of large amounts of Galactic diffuse gas invisible to the 
interferometer as in clouds G33.4 — 8.0 and G44.8 — 7.0, these values set an upper limit to 
the error possibly introduced by IMMERGE in the process of merging the VLA and GBT 
data. For most of the clouds this error is only a few per cent, indicating that the resulting 
cubes are scaled accurately. 



4.4. Converting to Galactic Coordinates 

The final image cubes were regrided to Galactic coordinates as the last step of the data 
reduction. To ensure a smooth transformation the pixel size was decreased. The final cubes 
are 512 x 512 with a pixel size of 5"7 over the 49' diameter field. The synthesized beam 
sizes and resulting gains for the data cubes are given in Table [3| 


4.5. The Noise Pattern 

The detection threshold for the 21-cm line emission depends on the noise and its 
distribution across the map areas. Examining the noise distribution is also a useful test 
of the effects of the merger of the GBT and VLA data. We have measured the noise in 
the cubes at various stages of the analysis. An essential step of the VLA data reduction 
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is the primary beam attenuation correction, PBCOR, which involves scaling the data to 
correct for the attenuation from the primary beam response of the 25m dishes of the VLA. 
This correction must be done before using IMMERGE to determine / ca i as described in 


section |4~T| The correction for the primary beam response produces images whose noise 
is a strong function of radius from the pointing center. This is illustrated in Fig. [6j which 
shows the VLA+GBT column density map for G26.9 — 6.3. The high noise near the edges 
of the map shows the effect of the primary beam correction on the interferometer data, 
which persists after merger with the GBT data. (The GBT maps have nearly uniform 
sensitivity everywhere.) The noise distribution has a characteristic shape a(r) as a function 
of radius, r, from the field center. The noise at the map edge is so high that it dominates 
the brightness scale of the image. 

To set a robust noise threshold for detection of the line, and to determine the errors 
in measured parameters based on the data, we need to understand the function er(r). The 


VLA primary beam gain factor is approximated in PBCOR by (Perley 2000): 


G(X) = 1 + G,X 2 + G 2 X 4 + G 3 X 6 , 
X = fr, 


( 2 ) 


where / is the frequency of observation in GHz and r is the angular radius in arcmin. The 
coefficients G{ are roughly constant for each VLA band; for our observations their values 
are: 


G x = -1.343 x 10" 3 , 

G 2 = +6.579 x 10“ 7 , (3) 

G 3 = -1.186 x 10“ 10 . 

PBCOR divides the spectrum at each point by G(fr ) from eq. thus amplifying both 
signal and noise, since 0 < G(fr ) < 1. The various clouds have slightly different values 
of the center frequency, /, depending on their radial velocities, and different values of (To, 
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which are almost entirely due to the noise in the VLA data, as explained in the next section. 
For setting detection thresholds and computing errors in the column density and mass at 
different points in each cloud, we use: 

cr(r) = <j 0 • [1 + Gi(fr) 2 + G 2 (fr) 4 + G 3 (fr) 6 } (4) 

Figure [7] shows the measured values of <r(r) using channels with no line emission as a 
function of distance from the center of the map shown in Fig. [6} The curve is the prediction 
of eq. [4], showing the effect of the primary beam correction applied to the VLA data, as 
it appears after merging with the GBT data. It is clear that noise from the VLA data 
completely dominates the noise in the final cubes, as the points in Fig. [7] are well described 
by the curve. For each cube we measure empirically the noise level at the field center, cr 0 , 
using off-line channels, and construct a noise function cr(r) similar to that shown in Fig. [7| 
These are given in Table [I] in Kelvins and the equivalent error in AVt for a 25 km s -1 wide 
velocity interval and the channel width of 0.64 km s -1 . 


4.6. Noise Amplitude 

Because of the angular resolution difference, noise from the GBT has little influence on 
the noise in the final cube, but rather appears as a systematic error in flux measurement 
over areas of a size comparable to the GBT beam. Rms noise values in the GBT cubes are 
0.08 - 0.14 K, a factor 3-4 times smaller than the noise at the very center of any of the final 
H I maps. Thus the dominant noise in the final data comes from the VLA, and errors due 
to the GBT noise can be neglected. 

We have processed the data in a number of non-trivial ways so it is important to 
check that the final noise level is reasonable, and is consistent with the noise in the VLA 
data at earlier reduction stages. Examination of the noise in line-free channels in the dirty 
continuum-subtracted cubes and in the clean cubes before PBCOR for four clouds is given 
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in Table [5| Comparing (cr c iean) for each cloud from Table [5] with a 0 in Table [5] we conclude 
that despite a significant number of processing steps following the cleaning of the VLA 
image, the rms noise value remains virtually unchanged, except by the correction for the 
main beam gain of eq. 


Two conclusions can be drawn: 1) because rms noise values are good indicators of the 
finest scale of the image, the fact that they do not change much from the VLA to the final 
data cube shows that our procedure of recovering the short spacings from the GBT data 
has not distorted the small-scale structure of the interferometric data; 2) the values of <To 
measured in § 4.5 and the noise pattern of eq. (jdj) are indeed valid indicators of the rms 
noise in the final data and can be used with confidence. 


4.7. A Noise Threshold 

Using eq. Q we can establish a noise threshold for every pixel in a cube. If there is no 
emission > 3cr(r) over the channel range of interest the spectrum is flagged and not used 
for making column density maps or other types of analysis. Figure [8] shows the same data 
as Figure |6j only with pixels blanked below the 3cr(r) level, leaving only emission detected 
significantly above the noise. 


5. Results 


5.1. Column Density Maps, Mass Profiles, Spectra 


Fig.|9| shows in comparison both full and thresholded column density maps of 
G16.0 + 3.0 based separately on the GBT and VLA data, and final VLA+GBT results. 


Fig. [10] gives the corresponding mass profiles. Finally, Fig. [TT| presents a summary of 
G16.0 + 3.0: thresholded VLA+GBT maps in relation to full GBT images, mass profile of 
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the VLA+GBT image, and spectra toward the Nhi peak of the clouds. Caption to the 
latter figure contains comments about the structure of the cloud. 


Figures [12- 38 present the same for the remaining 9 clouds. 


Measured cloud properties are summarized in Table [6} The values of line brightness 
temperature, velocity and line width are derived from the Gaussian decomposition, sampled 
toward the position of the peak Nhi • Errors are la from the Gaussian fit. Four of the 
clouds have spectral lines that require two Gaussians. The values of Nhi in col. 7 are 
integrals over the relevant velocity range (given in caption to each Figure) and are always 
close to the value from the sum of the Gaussian fits. We do not know the optical depth and 
so we cannot correct for self-absorption in the 21cm line, therefore the column densities and 
masses calculated in this paper are all lower limits. 


The mass profiles were constructed in the following fashion. The mass was sampled 
over a set of annuli of varying radius increment but equal area, all centered at the main 
column density peak of the map. With equal areas, the mass in each annulus is proportional 
to the average Nhi with the same proportionality coefficient and thus the points can be 
plotted simultaneously with two sets of axes: distance-dependent mass M vs. linear radius 


and distance-independent ( A H i ) vs. angular radius (see Fig. 10 and its counterparts, upper 
panels). The vertical error bars show the cumulative error due to nois^J the horizontal ones 
show the average beam radius rb ea m = ^(B ma jB min ) 2 , with B ma j and B rmn for the VLA 
given in column 2 of Table [3} The GBT values are the same for all clouds: 10.2'x 9.7', 
slightly increased compared to the original GBT beam due to pre-IMMERGE processing. 


Since there is a significant freedom of choice for such annulus sets, we have selected a 
fixed maximum radius r n of the sequence, the same for all clouds and equal to 0?319, which 


2 Both the mass and the linear radius measurements are also subject to the distance 
uncertainty (see Table |Tj) not shown in these plots. 
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covers most of the VLA primary beam, with the exception of a small outer portion having 
the highest noise. Then we required the width of the largest annulus to be equal to r'| )Cam , 
i.e. at the resolution of the map. Designating n the number of annuli and rq the smallest 


radius of the set, we arrive at the following equations: 

Ui f'n—l ^beam) (5) 

r i — r i-i = r\, ie[2,n], ( 6 ) 

rl-rU = ^ i G [2, n], (7) 

r i=irj, i G [2, n], (8) 

Solving eq. [5] together with eq. [8] at i = n and i = n — 1 we get: 

D [^’beam(2?' n ^beam)]^ • (9) 


Using eqns. [8] and [9] we have constructed the appropriate set of annuli for each cloud. 

Comparing line parameters of the GBT-only observations of Table [l] with the high 
resolution GBT+VLA results in Table [6] we find identical mean velocities, with a difference 
of 0 ± 2 km s _1 . As expected given the small angular structure revealed in the maps, 
the lines in the final maps are brighter by factors that range from 1.6 to 11 (for the very 
compact G16.0 + 3.0), with a median value around 4. These ratios are smaller than might 
occur: hydrogen clouds with sizes < V should appear ~ 100 times brighter to the VLA 
than the GBT, so it seems that the major structures in the clouds have been resolved in 
the current data. The line widths have a much smaller variation with the increased angular 
resolution. For lines that appear narrow to the GBT it is typical that they are even narrower 
in the combined data by around 20%, suggesting that the higher resolution observations 
are revealing colder or less turbulent material. But the broadest lines as measured with the 
GBT are sometimes even broader in the combined data, clouds G24.3 — 5.3 and G25.2 + 4.5 
being examples. 
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5.2. Estimating Masses, Sizes and Densities of the Clouds 


Table [7] gives the derived properties of the clouds for their adopted distances. Errors 
on derived parameters are dominated by the the distance uncertainties, whose estimates are 


discussed in £2.1 


In most cases the clouds do not display clear boundaries. In order to determine 
meaningful sizes and masses we have employed contours of constant iV H i. In each case a 
contour was chosen to encompass most of the visible cloud structures in the VLA+GBT 
maps (bottom left panels of Figures [9| [l2]etc.). The resulting contours are shown in the left 
panels of Figures [39| — 142} The H I masses inside these contours for each cloud are listed in 
column 4. For comparison we have also drawn contours in the GBT maps at much lower 


A^h i values. The resulting contours are shown in the right panels of Figures [39] - [42] and the 
derived properties are listed in Table [8} 


For each contour the major axis (the longest distance between two contour points) and 
the minor axis (the longest distance between contour points in the direction, perpendicular 
to the minor axis) are determined. The length of these axes are listed as D ma j x D min in 
column 5 of Table [7| and column 3 of Table (HI 


The volume density in column 6 of Table [7] and column 4 of Table [8] is the total H I mass 
divided by the volume given by where the effective radius r e = \{D ma jD m i n )^. 


5.3. Cloud Cores 

Table [9] gives estimates of the properties of the cloud cores, the denser regions of each 
cloud, determined by analyzing the H I emission at a higher value of IVhi around the 
column density peak of each cloud (column 2). The mass, size and number density of the 
enclosed area is given in cols. 3-5. The FWHM in Table [6] can be used to limit the kinetic 
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temperature — this value, Ti imit is given in col. 6. Finally, multiplying values in columns 5 
and 6 we can get a rough estimate of the core pressure in each cloud (column 7). For the 
clouds with two velocity components the total number density was split proportionally to 
the corresponding column densities. 


5.4. Discussion 


We have produced high angular-resolution 21cm H I maps of ten clouds that he in 
the boundary between the disk and the halo in the inner Galaxy. This paper presents the 
data and the reduction methods necessary to insure accurate results. We defer a detailed 
discussion of the cloud properties to a separate publication, but some general comments 
can be made, for these clouds are unique samples of the neutral interstellar medium. 

The disk-halo clouds, with masses of many hundreds of M© and locations many 
hundreds of pc from the Galactic plane, are orders of magnitude denser than their 
surroundings. 

The following rough estimates are made with the assumption of a spherical cloud of 
pure monoatomic hydrogen at constant number density n, its mass density p = By 

the cloud “Size” we understand its diameter and take that T iimi t. determined based on the 
emission line FWHM, represents its true kinetic temperature T. 


Under these conditions the sound-crossing time can be expressed as: 

Size „ ,, Size 


t 


sound 


2 Myr 


( 10 ) 


c, FWHM 

For clouds’ dense cores of sizes ~ 10 pc this is just a few Myr. But for whole clouds, with 
sizes ~ 100 pc, t S ound may reach 60 Myr. 


On the other hand, the free-fall time of gravitational collapse: 

tfr — \l — — ~ 50 Myr ■ n~ U 


32 Gp 


( 11 ) 
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For the highest cloud core density (n) ~ 10 cm -3 , tg ~ 16 Myr, but for other cores it 
is 20-50 Myr and for whole clouds this time scale often exceeds 100 Myr. So in all cases 
t sound <lff and thus the clouds are not gravitationally bound. 


It is instructive also to compare t SOU nd with the time t z of vertical fall to the Galactic 
plane at cloud’s location. Using a simple analytical expression of Wolfire et al. (1995b) for 
the ^-component of the Galactic gravitational acceleration g z (R , z), one can see that for 
\z\ < 200 pc the acceleration is close to linear: g z (R, z ) ~ g z (R ) • -2. The values of g' z range 
from (10.8 Myr) -2 at R = 8.5 kpc to (5.4 Myr) -2 at R = 2.5 kpc. Since motion with such 
acceleration is harmonically periodic, the free-fall time is just | g' z -l,/2 . For higher \z\ where 
g z (R, z ) flattens, these values can be used as lower limits: 


7r 


t* > TT fe'(fl)) 


( 12 ) 


equal to 8 - 17 Myr for 2.5 kpc < R < 8.5 kpc. In Pidopryhora (2006) a more precise 
ballistic calculation was done for 2 kpc < R < 5.3 kpc and a much greater distance from 


the plane z ~ 3.4 kpc using Walter Dehnen’s GalPot package (Dehnen & Binney 1998), 
obtaining t z ~ 30 Myr, which can be used as an upper limit for all the disk-halo clouds. 


For the cloud cores the condition t SOU nd W t z is true, so the clouds have time to respond 
to local physical conditions. The clouds must have come to internal pressure equilibrium, 
although the pressure may have contributions from turbulence on a range of scales, and 
possibly magnetic field and cosmic ray pressures as well. But since t sound > t z for the larger 
cloud structure, the overall density distribution of the cloud does not have time to dissipate 
in the low pressure of the halo over the time it takes the cloud to rise or to fall back to the 
disk. 


Similarly one can estimate the clouds’ Jeans masses as 

3 1 

( ^knT\ 2 / 3 \ 2 i 

—— -— ) « 9400 M q ■ FWHM 3 n - U (13) 

Gm H J V 47i 7V 

which is two to three orders of magnitude larger than their observed gas masses. 
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All this implies that the clouds are dynamic entities whose properties must reflect 


their history as well as conditions at their current locations (e.g., Koyama & Ostriker 


2009; Saury et al. 2014). While the appearance of many of the disk-halo clouds suggests 


that they are interacting with their local environment producing the steep gradients in 
Ah i or asymmetric shapes, we find little correlation between location in the Galaxy and 
fundamental cloud properties, with one clear exception, G16.0 + 3.0 discussed below. 

In theories of the interstellar medium the local pressure is often the controlling factor 
in the structure of neutral clouds, and over a wide range of conditions in the Galaxy 
it is expected that neutral clouds could consist of two phases, one warm and one cold 


(Field et al. 1969; McKee & Ostriker 1977; Wolfire et al. 1995a; Jenkins 2012). Four of 


the clouds studied here have line profiles indicating the presence of two components with 
different temperatures at the location of their peak Ah i . It is unlikely that this results from 
confusion of unrelated material as all two-component clouds are located at \z\ > 590 pc 
from the Galactic plane. The components typically have similar but not identical values of 
Vlsr- The narrower line component contains between 30% and 50% of the total A H i, about 
the mass fraction expected from some simulations (Saury et al. 2014). However, clouds may 
have two phases not easily separable in their emission profiles as the cold gas may have 
large velocity fluctuations that blend it with the warmer emission (Vazquez-Semadeni 2012 


Saury et al. 2014). 


Cloud G16.0 + 3.0 is particularly interesting, as it may he within the area around the 


Galactic nucleus excavated by a hot wind: the “Fermi Bubble” (e.g., Bland-Hawthorn & 
Cohen 2003; Su et al. 2010). The boundaries of the region effected by the wind are not well- 
defined, especially at low latitudes, but the G16.0 + 3.0 cloud has such distinctive properties 
as to suggest that its environment is different from that of the other clouds. Within the hot 


wind the pressure is many times larger than the typical ISM (Bland-Hawthorn & Cohen 


2003; 

Carretti et al. 2013 

); this is expected to force clouds into a purely cold phase ( 

Field et 
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al. [1969; Gatto et al. 2015). The properties of G16.0 + 3.0, with its high density, narrow line 
width, and compact structure, are consistent with this interpretation. Moreover, G16.0 + 3.0 
is notably smaller then the other clouds, but its size is similar to that of the H I clouds 


found to be entrained in the nuclear wind (McClure-Griffiths et al. 2013). In contrast, the 


nearby cloud G17.5 + 2.2 is nearly indistinguishable from the other clouds studied here. It 
might be as close as 135 pc to G16.0 + 3.0, but is more likely at least 1 kpc away given 
the uncertainties in our assignment of distances. Neither cloud shows kinematic anomalies 


suggesting that they have been accelerated by the nuclear wind (McClure-Griffiths et al. 


2013) 


We can also compare the properties of our disk-halo clouds with the disk clouds of 


Stil et al. (2006) observed at similar resolution, but all at z ~ 0. Only one of their clouds, 


59.67 — 0.39 + 60 has an analogous morphology and size, with a few times larger mass and 
average density compared to the clouds of our survey. The properties of all their other 
clouds are close to our cloud cores, only in some cases displaying a few times larger average 
density. 

These disk-halo H I structures allow us to study unconfused interstellar clouds in a 
variety of locations; this may lead to a better understanding of physical conditions that have 


hitherto been manifest only in ensemble averages. When the GASKAP survey (Dickey et 


al. 2013) is done, and even more when the full SKA is ready, studies like this one will reveal 
many more such clouds, and the techniques developed here will be useful for understanding 
their properties. 


This research was supported in part by the Australian Research Council through grant 
DP110104101 to the University of Tasmania. YP has started working on this research 
project while he held a predoctoral fellowship at the NRAO and then did a part of it while 
being employed at the Joint Institute for VLBI in Europe (Dwingeloo, the Netherlands). 
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Table 1. Cloud Properties from GBT Observations 


Name 

dV 

Vlsr 1 

FWHM 

Distance 

z 


(K) 

(km s -1 ) 

(km s -1 ) 

(kpc) 

(P c ) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 


G16.0 + 3.0 

1.1 

+143 

5.4 

8.2 ±1.0 

±430 ± 50 

G17.5 + 2.2 

1.6 

+139 

16.2 

8.1 ±1.1 

+310 ± 40 

G19.5- 

3.6 

3.5 

+121 

13.1 

8.0 ±1.6 

-500 ± 100 

G22.8 + 4.3 

3.1 

+137 

7.8 

7.9 ± 1.4 

+590 ± 100 

G24.3 - 

5.3 

1.3 

+124 

13.5 

7.8 ±1.5 

-720 ± 140 

G24.7 - 

5.7 

1.8 

+127 

8.9 

7.8 ±1.5 

-770 ± 150 

G25.2 + 4.5 

2.6 

+147 

10.4 

7.7 ± 1.6 

+610 + 130 

G26.9 - 

6.3 

3.4 

+123 

8.8 

7.6 ±1.7 

-840 ± 190 

G33.4 - 

8.0 

1.1 

+102 

5.5 

7.2 ±2.0 

-1000 + 280 

G44.8 - 

7.0 

3.3 

+94 

9.0 

6.1 ±2.5 

-740 ± 300 
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Table 2. Parameters of the VLA Observations 


Cloud 

Pointing 

Coordinates 

(J2000) 

Integration 

Time, 

minutes 

Central 

Velocity, 

km s -1 

Calibrators 

G16.0 + 3.0 

18 /l 08 m 42.850 s -13°35T7" 

131 

+140 

1833-210, 

1331+305 

G17.5 + 2.2 

18 /l 14 m 46.672 s -12°39'31" 

105 

+135 

1833-210, 

1331+305 

G19.5 — 3.6 

18 ft 39 m 46.681 s -13 0 32'41" 

90 

+120 

1833-210, 

1331+305 

G22.8 + 4.3 

18 /l 17 m 20.803 s -06°59'31" 

143 

+120 

1812-068, 

1331+305 

G24.3 - 5.3 

18 /l 54 m 52.557 s -10°03'07" 

222 

+125 

1833-210, 

1331+305 

G24.7- 5.7 

18 /l 57 m 05.347 s -09°52'41" 

64 

+125 

1833-210, 

1331+305 

G25.2 + 4.5 

18 /l 21 m 24.210 s -04°49'06" 

172 

+145 

1812-068, 

1833-210, 





1331+305 

G26.9 - 6.3 

19 /l 03 m 07.240 s —08° 17'17" 

215 

+122 

1833-210, 

0137+331, 





1331+305 

G33.4 — 8.0 

19 ft 20 m 52.118 s —03°16'43" 

102 

+102 

1939-100, 

1331+305 

G44.8 - 7.0 

19 fe 38 m 25.318 s +07°17T6" 

105 

+94 

1941-154, 

1331+305 
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Table 3. Synthesized beam sizes 


Cloud 

FWHM 

(") 

Gain 

(K/Jy) 

(FWHM) 

(P c ) 

G16.0 + 3.0 

89.5 x 52.6 

129 

2.7 

G17.5 + 2.2 

85.7 x 56.5 

126 

2.7 

G19.5 — 3.6 

84.6 x 58.0 

124 

2.7 

G22.8 + 4.3 

72.9 x 56.1 

149 

2.4 

G24.3 - 5.3 

75.6 x 56.2 

143 

2.5 

G24.7- 5.7 

77.0 x 59.8 

132 

2.6 

G25.2 + 4.5 

71.2 x 55.6 

154 

2.4 

G26.9 - 6.3 

73.5 x 57.3 

145 

2.4 

G33.4 - 8.0 

73.0 x 58.3 

143 

2.3 

G44.8 - 7.0 

68.5 x 57.7 

154 

1.9 
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Table 4. Noise levels at the field center 


Cloud 

^0 

(K) 

^(Ahr/) 1 
(10 18 cm -2 ) 

G16.0 + 3.0 

0.45 

3.2 

G17.5 + 2.2 

0.58 

4.1 

G19.5 — 3.6 

0.51 

3.7 

G22.8 + 4.3 

0.55 

3.9 

G24.3 - 5.3 

0.46 

3.3 

G24.7- 5.7 

0.48 

3.5 

G25.2 + 4.5 

0.42 

3.0 

G26.9 - 6.3 

0.42 

3.0 

G33.4 - 8.0 

0.43 

3.1 

G44.8 - 7.0 

0.41 

3.0 


1 In assumption of 40 channels of the total width of 25 km s 1 
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Table 5. Propagation of Noise Through the Imaging 


Cloud 

(C’dirty) 

(mJy/beam) 

(^dirty) 

(K) 

('^'cleari) 

(mJy/beam) 

(^clean) 

(I<) 

G16.0 + 3.0 

3.3 

0.42 

3.6 

0.46 

G26.9 - 6.3 

2.7 

0.39 

2.9 

0.42 

G33.4 - 8.0 

2.8 

0.40 

3.0 

0.42 

G44.8 - 7.0 

2.6 

0.41 

2.7 

0.42 
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Table 6. Measured Cloud Properties 


Cloud 

(i) 

R 

(kpc) 

(2) 

z 

(pc) 

(3) 

t 6 

(K) 

(4) 

Plsr 

(km s _1 ) 

(5) 

FWHM 

(km s -1 ) 

(6) 

Peak N hi 

(10 20 cm -2 ) 

(7) 

G16.0 + 3.0 

2.34 

+430 

12.3 + 0.3 

143.6 + 0.1 

4.7 ±0.1 

1.2 

G17.5 + 2.2 

2.56 

+310 

6.2 ±0.3 

138.3 ± 0.4 

18.2 + 1.0 

2.1 

G19.5 — 3.6 

2.84 

-500 

14.6 + 0.4 

118.8 + 0.1 

9.9 ±0.3 

2.9 

G22.8 + 4.3 

3.29 

+590 

4.9 ±0.4 

135.2 + 0.2 

5.8 ±0.6 

1.7 




3.2 ±0.4 

131.7 + 0.8 

24.3 + 2.2 


G24.3 - 5.3 

3.50 

-720 

2.7 ±0.1 

125.6 + 0.5 

17.4 + 1.1 

0.9 

G24.7 — 5.7 

3.55 

-770 

7.8 ±0.3 

128.0 + 0.1 

5.5 ±0.2 

1.0 

G25.2 + 4.5 

3.62 

+610 

6.3 ±0.2 

144.5 ± 0.2 

13.4 + 0.5 

1.6 

G26.9 — 6.3 

3.85 

-840 

13.0 + 0.4 

120.3 + 0.1 

3.0 ±0.1 

2.2 




5.2 ±0.3 

124.1 + 0.3 

14.5 + 0.6 


G33.4 - 8.0 

4.68 

-1000 

3.5 ±0.3 

101.9 + 0.2 

4.3 ±0.5 

0.8 




1.7 ±0.3 

101.0 + 0.8 

20.2 + 3.1 


G44.8 - 7.0 

5.99 

-740 

6.8 ±0.4 

93.4 + 0.1 

6.1 ±0.4 

1.7 




2.3 ±0.3 

99.4 + 1.1 

18.7 + 1.5 
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Table 7. Derived Cloud Properties 


Cloud 

R 

z 

M m 1 

Size 1 2 

(n) 


(kpc) 

(P c ) 

(M e ) 

(pc) 

(cm -3 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


G16.0 + 

3.0 

2.34 ±0.20 

±430 ± 50 

45i;s 

16 

X 

10 

(0-12) 

1.73 ±0.20 

G17.5 + 

2.2 

2.56 ±0.23 

±310 ± 40 

440:™ 

43 

X 

23 

(0-14) 

1.10 ±0.15 

G19.5 — 

3.6 

2.84 ± 0.42 

-500 ± 100 

83011™ 

67 

X 

28 

(0.20) 

0.80 ±0.16 

G22.8 + 

4.3 

3.29 ±0.29 

±590 ± 100 

550™ 

50 

X 

35 

(0.18) 

0.59 ±0.11 

G24.3 - 

5.3 

3.50 ±0.31 

-720 ± 140 

300l}ll 

59 

X 

29 

(0.19) 

0.33 ±0.06 

G24.7 - 

5.7 

3.55 ±0.30 

-770 ± 150 

290 +12 ° 

^au_ioo 

60 

X 

38 

(0.19) 

0.21 ±0.04 

G25.2 + 

4.5 

3.62 ±0.34 

±610 ± 130 

69011° 

58 

X 

48 

(0-21) 

0.37 ±0.08 

G26.9 - 

6.3 

3.85 ±0.36 

-840 ± 190 

60011° 

66 

X 

32 

(0.22) 

0.48 ±0.11 

G33.4 - 

8.0 

4.68 ±0.41 

-1000 ± 280 

160lso° 

38 

X 

28 

(0.28) 

0.36 ±0.10 

G44.8 - 

7.0 

5.99 ±0.50 

-740 ± 300 

84011° 

62 

X 

41 

(0.41) 

0.51 ±0.21 


Clouds’ sizes and masses determined based on contours shown in the left panels of 


Figs. G9-42 


fractional uncertainties in both dimensions are identical and are given by the value in 
parentheses. 
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Table 8. Wide-Field Cloud Properties Based on GBT Data 


Cloud 

m hi 1 

Size 1 2 

(n) 


(M & ) 

(P c ) 

(cm -3 ) 

(1) 

(2) 

(3) 

(4) 


G16.0 + 3.0 

83l?J 

47 x 43 (0.12) 

0.07 ±0.01 

G17.5 + 2.2 

750l 2 go 

57 x 39 (0.14) 

0.56 ±0.08 

G19.5 — 3.6 

1920lggo 

100 x 71 (0.20) 

0.25 ±0.05 

G22.8 + 4.3 

I 370 I 440 

103 x 70 (0.18) 

0.17 ±0.03 

G24 - 5 3 

II 5 OI 4 QQ 

158 x 59 (0.19) 

0.10 ±0.02 

G25.2 + 4.5 

1690lg3Q 

119 x 75 (0.21) 

0.16 ±0.03 

G26.9 - 6.3 

1570lg2Q 

151 x 85 (0.22) 

0.08 ±0.02 

G33.4 - 8.0 

540^260 

115 x 58 (0.28) 

0.08 ±0.02 

G44.8 - 7.0 

2560t4g 3 Q 

122 x 85 (0.41) 

0.19 ±0.08 


1 Clouds’ sizes and masses determined based on contours shown in the right panels of 


Figs. [39-42 


fractional uncertainties in both dimensions are identical and are given by the value in 
parentheses. 


3 Total for the group of three clouds, containing G24.3 — 5.3 and G24.7 — 5.7. 
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Table 9. Cloud Core Properties 


Cloud 

contour Nm 


Size 1 

K) 

Tlimit 

Tlimit ' (n c ) 


(10 20 cm -2 ) 

(M e ) 

(P c ) 

(cm -3 ) 

(K) 

(10 3 K cm -3 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 


G16.0 + 3.0 

0.2 

45 +12 

40 -10 

16 x 10 (0.12) 

1.7 ±0.2 

480 ± 20 

0.8 ±0.14 

G17.5 + 2.2 

1.5 

66i|? 

14 x 6 (0.14) 

6.7 ±0.9 

7200 ± 790 

48 ± 12 

G19.5 — 3.6 

2.2 

60ill 

10 x 6 (0.20) 

10.0 ±2.0 

2100 ± 130 

21 ±6 

G22.8 + 4.3 

1.2 

57±?1 

11 x 9 (0.18) 

4.5 ±0.8 







1.2 ±0.2 

730 ± 150 

0.9 ±0.3 





3.3 ±0.6 

13000 ± 2300 

42 ± 15 

G23.1 + 4.3 2 

0.9 

38±i® 

10 x 6 (0.18) 

6.4 ± 1.2 







2.2 ±0.4 

590 ± 200 

1.3 ±0.7 





4.2 ±0.8 

7300 ± 2200 

30 ± 14 

G24.3 - 5.3 

0.6 

45111 

15 x 14 (0.19) 

1.2 ±0.2 

6600 ± 830 

7.6 ±2.4 

G24.7- 5.7 

0.5 

£Q+22 

°°-18 

21 x 7 (0.19) 

2.3 ±0.4 

660 ± 50 

1.5 ±0.4 

G25.2 + 4.5 

1.2 

57 I 21 

15 x 7 (0.21) 

4.1 ±0.9 

3900 ± 290 

16 ±5 

G26.9 - 6.3 

1.3 

541 ll 

10 x 7 (0.22) 

7.2 ± 1.6 







2.5 ±0.5 

200 ± 13 

0.5 ±0.14 





4.7 ± 1.1 

4600 ± 380 

22 ±7 

G33.4 — 8.0 

0.5 

47l 3 ° 

24 x 8 (0.28) 

1.4 ±0.4 







0.4 ±0.1 

400 ± 90 

0.2 ±0.1 





1.0 ±0.3 

8900 ± 2700 

9 ± 5 
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Table 9—Continued 


Cloud 

contour A7n 

M H i 

Size 1 

( n c) 

Tlimit. 

Tii m it * (n c ) 


(10 20 cm -2 ) 

(M 0 ) 

(P c ) 

(cm -3 ) 

(K) 

(10 3 K cm -: 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

G44.8 - 7.0 

1.2 

54l| 

15 x 7 (0.41) 

3.9 ± 1.6 







1.9 ±0.8 

810 ± no 

1.6 ±0.8 





2.0 ±0.8 

7600 ± 1200 

15 ±9 


1 Fractional uncertainties in both dimensions are identical and are given by the value in 
parentheses. 

2 A small separate cloud at the high-longitude side of the G22.8 + 4.3 held. First velocity 
component: T& = 6.4 ± 0.9 K, Vlsr = 129.3 ± 0.3 km s -1 , FWHM = 5.2 ± 0.9; second 
velocity component: T& = 3.4 ± 0.6 K, Vlsr = 135.1 ± 1.7 km s -1 , FWHM = 18.3 ± 2.7. 
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Fig. 1.— Velocity-latitude image of H I emission from GBT observations after Lockman 


& Pidopryhora (2005). The cloud at 24.3 — 5.3 is marked with an arrow. The expected 


maximum Vlsr from Galactic rotation (Clemens 1985) is marked by the slightly curved 
vertical line. The location of G24.3 — 5.3 slightly beyond the maximum velocity indicates 
that it must he near the tangent point in its direction. Note that it is separated in position 
and velocity from other H I emission and is thus relatively unconfused. 


34 




























V (10 X) 



Fig. 2.— un-coverage of the VLA observations of one of the clouds, G26.9 —6.3. Two epochs 
of June 29 and July 24 2004 combined, total integration time 215 min., source declination 
- 8 °. 
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High Resolution Data (Jy) 


Fig. 3.— Linear fit to pixel-by-pixel comparison of single-dish and interferometer data of 
cloud G17.5 + 2.2 generated by Miriad’s task IMMERGE in the process of solving for the 
optimal / ca i value. The left panel is without tapering, the right panel uses tapering. Data are 
used from the position of the brightest line averaged over a region 1513 x 1513. The overlapping 
frequency range in the Fourier domain was set to be 20 - 50 meters « 95 - 240A. The values 
for the high resolution data were calculated from the interferometric data convolved with 
the single-dish beam to compensate for the difference in resolution. The values for the low 
resolution, single-dish data were scaled by / ca i . The value of / ca i determined from the left 
panel is 0.68, from the right is 0.62. In the tapered case it is clear that the assumption 
of linearity fails, while in the non-tapered case the bulk of the points close to zero favor a 
somewhat different slope. 
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Fig. 4.— Page 1 of 3 — Mean line profiles in brightness temperature T*,, at the center of each 
field for clean VLA data (orange), GBT data (purple) and results of running IMMERGE on 
these two images with / ca ,i =0.87 and no tapering (green). This sum is calculated over all 
pixels in an 80 by 80 pixel box. The comparison shows how well the total single-dish flux is 
recovered by IMMERGE. 
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Fig. 4.— Page 2 of 3 — Mean line profiles in brightness temperature T*,, at the center of each 
field for clean VLA data (orange), GBT data (purple) and results of running IMMERGE on 
these two images with / ca ,i =0.87 and no tapering (green). This sum is calculated over all 
pixels in an 80 by 80 pixel box. The comparison shows how well the total single-dish flux is 
recovered by IMMERGE. 
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Fig. 4.— Page 3 of 3 — Mean line profiles in brightness temperature T*,, at the center of each 
field for clean VLA data (orange), GBT data (purple) and results of running IMMERGE on 
these two images with / ca i =0.87 and no tapering (green). This sum is calculated over all 
pixels in an 80 by 80 pixel box. The comparison shows how well the total single-dish flux is 
recovered by IMMERGE. 
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Fig. 5.— Quantitative measure of the difference between GBT and IMMERGE spectra 
defined by eq. [l] The small differences indicate that the data reduction procedure has 
restored the missing flux in the line. 
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Fig. 6.— An example of an H I column density map with no noise threshold applied, showing 
the distinct radial increase in noise pattern arising from the VLA primary beam correction. 
This map of cloud G26.9 — 6.3 was made from 51 spectral channels covering the velocity 
range 108.47 < Vlsr < 140.68 km s -1 . 
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Fig. 7.— The rms noise value OT B ( r ) f° r the final G26.9 — 6.3 VLA+GBT cube as a function 
of angular distance r from the center of the field of view. The values of scattered points in 
this plot were calculated from the same data cube as the map in Fig. [6j only in a range of 
41 line-free channels 150.34 < Flsr < 176.11 km s -1 . The magenta curve fit was obtained 
by multiplying each point by the corresponding beam gain factor from eq. [2| averaging, then 
multiplying the resulting mean value <jq — 0.42 K by the same gain factor (eq. [3]). 
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Fig. 8.— H I column density map of G26.9 — 6.3 displaying the same data as Fig. [6j where 
spectra with values below 3<r(r) were flagged and not plotted, as explained in § |4.7[ The 
noise <j(r) was calculated using eq. (jdj) with values of <To and / from Table |4j The figure 
shows data over the interval 1 x 10 19 < Nhi < 2.2 x 10 2 ° cm -2 , where the lower limit is 
three times the column density noise at the center of the field. 
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Fig. 9.— Comparison of H I column density maps of G16.0 — 3.0, integrated over 36 
spectral channels in the interval 132.3 < Vlsr < 154.8 km s -1 . 3 channels at 139.4 - 
140.6 km s -1 were flagged due to RFI and replaced with a linear interpolation. Upper row: 
GBT data regridded to match the VLA. Middle row: VLA data. Bottom row: The result of 
IMMERGE, VLA+GBT. Left column: full data. The column density range of the full VLA 
and VLA+GBT maps (left middle and left bottom correspondingly) of this and all similar 
figures is fixed at —5 x 10 19 < A r H i < 1.7 x 10 2 ° cm -2 . Right column: The same maps as 
in left column, only thresholded by 3<r(r) of the VLA noise level, similar to Fig. [8j All three 
are shown over the same interval of column densities, 9 x 10 18 < A+j < 1.2 x 10 2 ° cm -2 , 
where the lower limit is three times the column density noise at the center of the field. This 
shows the GBT data (upper right) compared to both VLA signal and noise values. In this 
particular case, alone among the ten clouds under consideration, the GBT signal is never 
greater than the VLA 3a(r) level. 
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Fig. 10.— Radial mass distribution for G16.0 + 3.0. Top: mass and average column density 
per annulus, with sequence designed to ensure equal areas in each bin, as explained in § M 
Bottom: integral mass within circles of each radius. Red - GBT data, corresponding to 
Fig. [9] top left (only two of the error bars are shown to avoid clutter). Green - VLA data, 
corresponding to Fig. [9]middle left. Blue - VLA+GBT data, corresponding to Fig. [9]bottom 
left. 
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Fig. 11.— Summary of G16.0 + 3.0. Top left: 3cr(r) thresholded VLA+GBT map, same 
as Fig. [9] bottom right. Top right: GBT column density map integrated over the same 
velocity range, with the box marking the extent of the VLA map. Lower left: radial mass 


distribution, same as blue points in Fig. [10] top. Lower right: final spectrum toward the 
position of the peak AG 1 with the Gaussian whose parameters are given in Table [6] This is 
a very little cloud with no structure, unlike any other cloud in the group. It is the smallest 
cloud by far, has one of the narrowest line widths implying a kinetic temperature < 500 K, 


and has the largest mean density by a factor ~ 5 (see §5.4). About one half of the cloud 
mass is confined to the dense and very compact core. It may be in the region of the Galaxy 


excavated by a hot wind from the Galactic nucleus as discussed in §5.4| This isolated cloud 
is similar to the dense cores seen in other clouds, e.g., G26.9 — 6.3. 
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Fig. 12.— H I column density maps for G17.5 + 2.2, integrated over 48 spectral channels in 
the interval 122.8 < Vlsr < 153.0 km s -1 , as described in the caption to Fig. [9] 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 





























r[°] 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



Fig. 13.— Radial mass profiles for G17.5 + 2.2 as described in the caption to Fig. [TO} 
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Fig. 14.— Summary of G17.5 + 2.2 as described in the caption to Fig. [TTJ A cloud with 
an almost circular core ~ 30 pc in diameter that contains two peaks of similar Ahi and 
about half of the cloud’s mass. Surrounding the core are irregular structures of low Nm, 
the densest of which has a cometary appearance stretching to lower longitude and latitude. 
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Fig. 15.— H I column density maps for G19.5 — 3.6, integrated over 56 spectral channels in 
the interval 107.8 < Vlsr < 143.2 km s -1 , as described in the caption to Fig. |9j 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 































r[°] 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



Fig. 16.— Radial mass profiles for G19.5 — 3.6 as described in the caption to Fig. [lOj 
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Fig. 17.— Summary of G19.5 — 3.6 as described in the caption to Fig. [TTJ One of the 
two largest clouds in the sample with a diameter > 60 pc and a mass around 800 M 0 in 
the central region and more than 1000 M 0 in its extended envelope. The cloud’s boundary 
is sharpest away from the Galactic plane and there are two “tails” of less bright emission 
stretching back toward the plane. 
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Fig. 18.— H I column density maps for G22.8 + 4.3, integrated over 46 spectral channels in 
the interval 125.2 < Vlsr < 154.1 km s -1 , as described in the caption to Fig. |9| 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 
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Fig. 19.— Radial mass profiles for G22.8 + 4.3 as described in the caption to Fig. [TO} 
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Fig. 20.— Summary of G22.8 + 4.3 as described in the caption to Fig. [TTJ This cloud consists 
of several fragments, one of which (G23.1 + 4.3) has a size of 20 pc and a mass of 70 M 0 , 
similar to the small isolated cloud G16.0 + 3.0. The line at the peak Nm consists of two 
components, one broad and one narrow. 
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Fig. 21.— H I column density maps for G24.3 — 5.3, integrated over 70 spectral channels in 


the interval 107.0 < Vlsr < 151.4 km s 1 


as 


described in the caption to Fig. |9j 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 
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Fig. 22.— Radial mass profiles for G24.3 — 5.3 as described in the caption to Fig. [lOj 
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Fig. 23.— Summary of G24.3—5.3 as described in the caption to Fig. [IT} This cloud, together 
with G24.7 — 5.7 and a smaller one that lies between them, form a chain of three clouds 
connected by diffuse H I emission. G24.3 — 5.3 has no clear core, but rather several column 
density peaks of similar values spread over the body of the cloud. Its radial mass profile 
shows the lack of central concentration. In contrast to G16.0 + 3.0, which is a core without 
an envelope, this cloud is an envelope without a core. With a FWHM of 17.4 ± 1.1 km s -1 , 
its line is one of the broadest in the sample. 
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Fig. 24.— H I column density maps for G24.7 — 5.7, integrated over 66 spectral channels in 
the interval 103.7 < Vlsr < 145.6 km s -1 , as described in the caption to Fig. |9j 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 
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Fig. 25.— Radial mass profiles for G24.7 — 5.7 as described in the caption to Fig. [lOj 
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Fig. 26.— Summary of G24.7 — 5.7 as described in the caption to Fig. [TTJ A relatively small 
cloud that together with G24.3 — 5.3 and a cloud between them forms a chain of three clouds 
connected spatially and kinematically. Unlike its companion G24.3 — 5.3, it has a clear core 
with a narrow line of FWHM=5.5 km s -1 , but still has a similar range of column densities. 
The edge of the chain of clouds is quite sharp at higher longitudes. 
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Fig. 27.— H I column density maps for G25.2 + 4.5, integrated over 64 spectral channels in 
the interval 125.7 < Vlsr < 166.3 km s -1 , as described in the caption to Fig. |9| 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 
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Fig. 28.— Radial mass profiles for G25.2 + 4.5 as described in the caption to Fig. [TO} 
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Fig. 29.— Summary of G25.2 + 4.5 as described in the caption to Fig. [TTJ The VLA data 
reveal the considerable complexity of this cloud. There seems to be no single core, but 
instead many dense regions. The broad line with FWHM=13.4 km s -1 may indicate the 
blending of several kinematically distinct components within the If05 angular resolution of 
the map, equivalent to 2.35 pc at the adopted distance of the cloud. The cloud has a relatively 
sharp edge to higher longitudes with diffuse material spreading out to lower longitudes. The 
relatively broad mass profile and the GBT data suggest that the diffuse components of the 
cloud may extend well beyond the VLA field of view. 
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Fig. 30.— H I column density maps for G26.9 — 6.3, integrated over 51 spectral channels in 
the interval 108.5 < Vlsr < 140.7 km s _1 , as described in the caption to Fig. |9j 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 
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Fig. 31.— Radial mass profiles for G26.9 — 6.3 as described in the caption to Fig. [lOj 
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Fig. 32.— Summary of G26.9 — 6.3 as described in the caption to Fig. [TT} A cloud with a 
sharp boundary toward high latitude and toward the plane, and diffuse material spreading 
out on the opposite side. More than half of the mass is within 20 pc of the brightest point, 
but the diffuse material is responsible for the broad mass profile. Its spectrum at peak A H i 
shows two components with broad and narrow FWHM of 14.5 and 3.0 km s -1 , suggestive of 
a cloud with a two-phase structure. The line widths of the two components place an upper 
limit on their temperatures of 4600 K and 200 K. 
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Fig. 33.— H I column density maps for G33.4 — 8.0, integrated over 37 spectral channels in 
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the interval 92.3 < Vlsr < 115.5 km s 1 


as 


described in the caption to Fig. [9} 
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[10 19 cm 2 ] V HI [10 19 cm 2 ] N R1 [10 19 cm : 
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Fig. 34.— Radial mass profiles for G33.4 — 8.0 as described in the caption to Fig. [lOj 
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Fig. 35.— Summary of G33.4 — 8.0 as described in the caption to Fig. E3 This small 
isolated cloud of ~ 160 M e in its dense part is the most distant from the Galactic plane 
of any in the sample and has a complex structure. It has many column density maxima 
spread along a central ridge on the edge closest to the Galactic plane. The spectrum at the 
highest column density point has two components, one broad and one narrow, suggesting a 
two-phase thermal structure. 
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Fig. 36.— H I column density maps for G44.8 — 7.0, integrated over 73 spectral channels in 
the interval 79.8 < Vlsr < 126.2 km s -1 , as described in the caption to Fig. [9} 
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[10 19 cm 2 ] iV HI [10 19 cm 2 ] N R1 [10 19 cm : 
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Fig. 37.— Radial mass profiles for G44.8 — 7.0 as described in the caption to Fig. [lOj 
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Fig. 38.— Summary of G44.8 — 7.0 as described in the caption to Fig. [TT} In the spectrum 
a channel at 74.0 km s -1 is flagged due to an RFI spike. This large, massive cloud extends 
beyond the primary beam of the VLA. This, and the very small cloud G16.0+3.0 are the only 
clouds whose shape is nearly circular, though the brightest gas is concentrated in numerous 
small clumps. The spectrum at the ATi peak is asymmetric with a broad wing to higher 
velocity. This cloud has among the largest excess Vlsr beyond the expected V t of any in the 
sample. 
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Fig. 39.— Contours for determination of size and mass of clouds based on VLA+GBT (left) 
and GBT only (right). Top: G16.0 + 3.0, contour levels iVni = 2.0 x 10 19 cm -2 (VLA+GBT) 
and 5.5 x 10 18 cm -2 (GBT); bottom: G17.5 + 2.2, contour levels V H i = 5.6 x 
10 19 cm -2 (VLA+GBT) and 4.4 x 10 19 cm -2 (GBT). The dashed lines designate the axes: 
major (the longest distance between two contour points) and minor (the longest distance 
between contour points in the direction, perpendicular to the minor axis). 
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Fig. 40.— Contours for determination of size and mass of clouds based on VLA+GBT (left) 

and GBT only (right). Top: G19.5 —3.6, contour levels Nm = 5.5 x 10 19 cm -2 (VLA+GBT) 

and 3.9 x fO 19 cm -2 (GBT); middle: G22.8 + 4.3, contour levels V H i = 4.1 x 

10 19 cm -2 (VLA+GBT) and 2.6 x 10 19 cm -2 (GBT); bottom: G25.2 + 4.5, contour levels 

N m = 4.3 x 10 19 cm” 2 (VLA+GBT) and 2.2 x 10 19 cm” 2 (GBT). 
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Fig. 41.— Contours for determination of size and mass of clouds based on VLA+GBT (left) 
and GBT only (right). Top left: G24.3 —5.3 (VLA+GBT), contour at V H i = 3.1 x 10 19 cm -2 ; 
bottom left: G24.7 — 5.7 (VLA+GBT), contour at Vri = 2.3 x 10 19 cm -2 ; right: Group of 
G24.3 — 5.3 and G24.7 — 5.7 (GBT), contour at V H i = 1.9 x 10 19 cm -2 . 


76 


[10 19 cm 













8.9 



Fig. 42.— Contours for determination of size and mass of a cloud based on 

VLA+GBT (left) and GBT only (right): Top: G26.9 — 6.3, contours at iVni = 5.3 x 

10 19 cm -2 (VLA+GBT) and 1.1 x 10 19 cm -2 (GBT); middle: G33.4 — 8.0, contours at 

Ah i = 2.5 x 10 19 cm” 2 (VLA+GBT) and 1.2 x 10 19 cm” 2 (GBT); bottom: G44.8 - 7.0, 

contours at AThi = 5.5 x 10 19 cm -2 (VLA+GBT) and 3.1 x 10 19 cm -2 (GBT). 
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